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We study the superfluid-Mott-insulator transition of ultracold bosonic atoms in a one-dimensional 
optical lattice with a double-well confining trap using the density-matrix renormalization group. 
At low density, the system behaves similarly as two separated ones inside harmonic traps. At 
high density, however, interesting features appear as the consequence of the quantum tunneling 
between the two wells and the competition between the "superfluid" and Mott regions. They 
are characterized by a rich step-plateau structure in the visibility and the satellite peaks in the 
momentum distribution function as a function of the on-site repulsion. These novel properties 
shed light on the understanding of the phase coherence between two coupled condensates and the 
off-diagonal correlations between the two wells. 
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I. INTRODUCTION 

The experimental realization of trapped ultracold 
bosonic gases on optical lattices has opened up the pos- 
sibility of observing various quantum phases (e.g., su- 
perfluid and Mott insulator) and studying the quantum 
phase transitions between them in a well-controlled man- 
ner. Without the confining trap, the system will undergo 
a superfluid-Mott-insulator phase transition at integer 
fillings. In the presence of a confining trap, the super- 
fluid and Mott-insulator regions may coexist^ due to 
the inhomogeneity induced by the trap. 

Numerical calculations have been carried out for the 
systems without^ and with a harmonic trap^^^, us- 
ing the quantum Monte Carlo and density-matrix renor- 
malization group2ii£ (DMRG) methods. Interesting new 
features appear in the system with a single harmonic con- 
fining trap, which are associated with the competition be- 
tween the "superfluid" and Mott insulator regimes and 
characterized by the "kinks"— in the visibility of interfer- 
ence fringes with the increase of the on-site repulsion. 

Recently a bosonic Joscphson junction composed of 
two weakly coupled Bose-Einstein condensates in a 
macroscopic double-well potential has been experimen- 
tally realized^. It raises an interesting question con- 
cerning the quantum tunnelling effect in a double-well 
potential system in which the ultracold atoms on opti- 
cal lattice are trapped. In such a system, the inter-atom 
interaction and the competition between the "superflu- 
idity" and Mott insulator phases are expected to play a 
crucial role in shaping the inter-well tunnelling. 

In this paper, we intend to clarify this issue by studying 
a one-dimensional Bose-Hubbard model in the presence 
of a double- well potential using the DMRG. At low fill- 
ing (p = 1/2), we find the system simply behaves as two 
separated ones, each of which is in a harmonic trap with- 
out significant coupling between them. At a higher fill- 
ing (p = 1), however, the bosons in the two wells become 



strongly correlated and a rich phenomenon related to the 
evolution of the "superfluid" and Mott insulator regions 
as a function of the on-site repulsion is obtained. A de- 
tailed discussion on how to detect these features based 
on the visibility, momentum distribution, hopping cor- 
relation function, and other physical quantities is given. 
These novel properties shed light on the understanding of 
the phase coherence and off-diagonal correlation between 
the coupled condensates in a double- well interacting Bosc 
system on optical lattices. 

A cold atomic Bose gas on a one-dimensional optical 
lattice in the presence of a double-well potential trap can 
be described by the following Bose-Hubbard Hamiltonian 
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is the double-well trap with quadratic Vti and quartic 
Vt2 coefficients, and L is the number of sites. The hop- 
ping integral t is set as the unit of energy t = 1 , and U is 
the on-site repulsion. fij and fi.; are the bosonic creation 

and annihilation operators, respectively, and hi = fife is 
the number operator. 

An important parameter in characterizing the phase 
coherence in a superfluid-Mott-insulator transition is the 
visibility of interference fringes defined b y 12 d 3 
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the momentum distribution function 
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In the presence of a confining trap, it is useful to intro- 
duce the local density of bosons n, =< hi > and the local 
particle fluctuation^ 

Ki = [(hf) -nf], (5) 

to measure the inhomogeneity induced by the confining 
trap. To further quantify the spatial correlations within 
the trap potential, we have also calculated the hopping 
correlation function 

g(hj) = (aVaj) (6) 

between two lattice sites i and j, the total interaction 
energy E v , 

E v = Uj^niini - l)/2, (7) 

i 

the kinetic energy \Ek\, 

E k = -tJ2 (a\a l+1 + h.c.) , (8) 

i 

and the trapping energy Et 

E t = £Vr(iK (9) 

i 

respectively. 

II. SUPERFLUID-MOTT-INSULATOR 
TRANSITION 

In the following we present our DMRG calculations 
for the Bose- Hubbard model JT]). Compared with a har- 
monic trapi, the present system shows a richer and more 
complex structure, due to the coupling between the two 
wells. When the density is low enough such that there is 
no particle in the middle of the system, the system splits 
into two weakly coupled harmonic ones. However, if the 
density is sufficiently high that particles in the two wells 
are not well separated, the interaction between the wells 
will play a crucial role, which may change the behavior 
of the system dramatically. 

In the following, we will consider the system with 80 
sites in a double-well trap potential with Vri — — 0.024i 
and Vt2 = 2 x 10~ 5 i. The calculations are carried out 
using the DMRG with open boundary conditions, for two 
different filling factors p = 1/2 and p = 1. 

A. Filling p = 1/2 

AT the filling p = 1/2, the system will not reach a 
Mott-insulator state, even in the Tonks-Girardeau (TG) 
limit U — > oo, without a trap potential Vr- In the pres- 
ence of the double- well potential, however, two separated 
Mott-insulator domains with the density n, = 1 appear 
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FIG. 1: Density profiles m and the corresponding particle 
fluctuation k, for p = 1/2. The profiles for U = 5.8t and 
U — 6.2t nearly coincide. 

inside the two wells at large U, due to the low density of 
the system. Fig. [1] shows the density profile and corre- 
sponding particle fluctuation at various values of U/t. At 
small U /t, two separate superfiuid regions appear in the 
two wells. When U is increased, the superfiuid regions 
shrink until the Mott-insulator domains emerge when U/t 
is larger than 6.8. 

To further quantify this transition, let us define the 
integrated density 

h 

N d = ^m (10) 

i—li 

in a range [li, 1%]. In the present system, we choose l\ = 
11 and I2 = 22 in the left-side well. As shown in Fig[2](a), 
Nd is saturated to a plateau of Nd = 12 corresponding to 
a local Mott-insulator state (rii = 1) as U/t is increased 
beyond 6.8. 

Fig[2ja) shows the visibility v. Two kinks can be 
observed. The first one (less evident) occurs around 
U/t = 5.5, while the second one around U/t = 6.4. As 
U/t is increased between 5.8 and 6.2, there is almost no 
change in the density profile. This is due to presence of 
the emergent Mott-insulator regions surrounding the su- 
perfiuid ones in the two wells. The atoms in the central 
superfiuid regions can transfer to the outer superfiuid re- 
gions at larger U/t. This will eventually exhaust the su- 
perfiuid regions in the two well centers. At U/t = 20, the 
flat Mott-insulator plateaus, with vanishing local parti- 
cle fluctuation, centered around the two wells are clearly 
seen. 

As shown in Fig. ® the visibility with the correspond- 
ing Et, Ejj and the ratio 7 = \Ejj/Ek\, exhibits a char- 
acteristic kink structure. The quick decrease of Ejj in 
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FIG. 2: (a) Visibility and integrated density Nd; (b) 
Interaction energy Eu and trapping energy Et\ (c) Ratio 
7 = \Eu/Ek\ of interaction to kinetic energy, as functions of 
U/t, for the filling factor p = 1/2. 



FIG. 4: (a) Visibility v and integrated density Nd from ii = 11 
to h — 30, (b) Kinetic \Ek\ and total energy Es, (c) Ratio 
7 = \Eu/Ek\ of interaction to kinetic energy, as functions of 
U/t for p=l. 
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FIG. 3: Momentum distribution function S(k) for the system 
shown in Fig. [1] 



tion function S(k) for the system shown in Figfl] When 
U/t is small, the system is in a supcrfluid state and a sin- 
gle narrow peak appears at zero momentum in S(k). For 
large U/t, S(k) is broadened. The quantitative change 
in S(k) is best reflected in the kink structure of the 
visibility v. In he Tonks-Girardeau limit, we find that 
Smax = 1.04 and v = 0.44. The non-zero visibility v re- 
veals the presence of the superfiuid regions surrounding 
the Mott-insulator plateau regions (Fig. [T|). 

Finally, we note that at p = 1/2, the distribution of 
bosons remains well separated in the two wells [Figfl]. 
Therefore, it can be approximately treated as two single 
wells with weak coupling. Indeed, the results presented 
above are consistent with those obtained in a harmonic 
trapi. 



the interval U/t = 6.2 — 6.8 shows a great suppression of 
the double occupancy in the two wells. It is a manifesta- 
tion of the transition from superfiuid to Mott-insulator 
in the two wells. This occurs even though the total en- 
ergy increases continuously with the increase of U ft. By 
further increasing U /t, the density profile remains almost 
unchanged. 

Fig. [3] shows the corresponding momentum distribu- 



B. Filling p = 1 

At p = 1, bosons inside the two wells cannot be sepa- 
rated for sufficiently large U . The system is in a Mott- 
insulator state without the trap in the TG limit. 

FigfJJa) shows the visibility v and the integrated den- 
sity N d {h = 11, l 2 = 30) defined in Eq.^. The evo- 
lution of these two quantities exhibit a rich step-plateau 
structure compared with the case of p = 1/2. When U is 
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FIG. 5: Density profiles m and the corresponding local par- 
ticle fluctuation m for p = 1. 



small (about < lOi), the central Mott-insulator domain 
has not formed. In this parameter regime, the plateaus 
in v and Nd are quite narrow, indicating that the Mott- 
insulator domains surrounding the superfluid ones are 
not robust and the bosons can escape from the superfluid 
regions in the two wells. This feature is well reflected in 
the density profiles and local particle fluctuation around 
the step between U/t = 9.9 and 10 as shown in Fig|5]Ja). 

When U/t is further increased, the central region of the 
system shows a flat Mott-insulator plateau with n, = 1 
[Fig. [5]. This leads to a sudden decrease in v and Nd- 
Due to the large on-site repulsion U and the existence of 
the central Mott-insulator regions, bosons can only tun- 
nel from the inner superfluid region to the outer ones. 
Whenever a boson escapes from the superfluid region, the 
visibility v as well as the integrated density Nd shows a 
sharp decrease. When U/t is larger than 22, the super- 
fluid region disappears compctcly (Fig|5Jb)). In the TG 
limit, the whole system is in a Mott-insulator state and 
Kj = 0. 

As shown in Fig. BJb), the abrupt changes are also 
present in \Ek\ at the steps of v, although the total en- 
ergy Es increases smoothly with U/t. The steps in \Ek\ 
can be understood as the abrupt redistribution of the 
boson density between the superfluid and Mott-insulator 
regions. During this process, the superfluid region be- 
comes smaller, while the Mott-insulator becomes larger. 
The monotonic decrease oi\Ek \ with increasing U/t is due 
to the widening of the Mott insulator regions in which the 
kinetic energy is suppressed. In the plateau regions, the 
increase in 7 is also due to the reduction of the kinetic 
energy. 

Fig. [5] shows the momentum distribution function 
S(k). S(k) exhibits different features at different U/t 
regimes. Furthermore, S(k) is correlated with v. When 
the satellite peaks appear in S(k), v undergoes a sharp 
increase; while v undergoes a sharp decrease when the 
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FIG. 6: Momentum distribution function S(k) at different 
values of U/t for p = 1 . 



satellite peaks disappear. 

More detailed study demonstrates that the satellite 
peak structure in S(k) is related to the existence of a 
strong off-diagonal long-range correlation between the 
two wells. This can be seen from Fig. [3 in which the 
hoping correlation function <?(16, i) between site 16 and 
site i is plotted. If the system has an off-diagonal long- 
range correlation, the satellite peaks are present in S(k). 
If the system does not have the off-diagonal long-range 
correlation, the satellite peaks disappear. These features 
are different from those in a harmonic trap£. All the 
above features appear before the formation of the central 
Mott-insulator plateau. In the presence of the central 
Mott-insulator plateau, however, the satellite structure 
of S(k) as well as the long-range off-diagonal correlation 
between the two wells disappears. This suggests that 
one can judge whether or not the off-diagonal long-range 
correlation exists from the measurement of momentum 
distribution functions. 



III. SUMMARY 

To summarize, we have explored the supcrfluid-Mott- 
insulator transition in a double- well trapped atomic Bose 
gas in a one-dimensional optical lattice using the DMRG. 
We find that this transition is well characterized by the 
visibility v, which exhibits a series of kink structures as- 
sociated with the redistribution of bosons between local 
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shows a number of characteristic features in the local in- 
tegrated density, the density profile, local compressibility, 
and the momentum distribution. 

At low density (e.g. p = f/2), the system behaves 
like two weakly coupled harmonic traps. With the in- 
crease of the density, the system begins to develop novel 
properties unique to the double well potential associated 
with a strong quantum tunneling effect between the two 
wells. At high density (e.g. p = 1), a series of steps 
and plateaus appear in the visibility and other physical 
quantities. These features are related to the abrupt redis- 
tribution of the bosons in the local superffuid and Mott- 
insulator regions. In addition, the presence of the satel- 
lite peaks in the momentum distribution function indi- 
cates that the strong off-diagonal long-range correlation 
between the supcrfluid regions are separated by the Mott 
insulator regions. Therefore, it is an effective probe for 
experiment to identify the off-diagonal long-range corre- 
lation and phase coherence due to the tunnelling between 
the two condensates on the optical lattice. 
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the on-site repulsion U . The evolution of the system also 
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